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Abstract 

We show that periodically folded graphene sheets with enhanced spin-orbit interaction due to 
curvature effects can carry spin polarized currents and have gaps in the electronic spectrum in the 
presence of weak magnetic fields. Our results indicate that such origami-like structures can be used 
efficiently in spintronic applications. 
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Often hailed as a wonder material due to its impressive physical properties [T] , graphene 
has opened several venues of basic science exploration and it is a material that has a tremen- 
dous technological potential. Despite its potential applicability, the lack of a bandgap is a 
well known limitation that currently prevents the use of graphene in digital electronic appli- 
cations [2]. Different strategies have been attempted to remedy this shortcoming, namely, 
by quantum confinement in the form of nanoribbons and quantum dots [31 H], by stacking 
graphene sheets in bilayers in the presence of a perpendicular electric field |5], |6], by strain- 
engineering its electronic structure [7HS] , or by simply chemically doping the graphene sheets 
[TUt [TT] . Unfortunately, these attempts have so far failed to produce technologically relevant 
semiconducting graphene due to several difficulties that go from the small size of the gaps 
they produce to the disorder that they introduce. 

On a different front, the field of spintronics appears as one of the most promising areas 
for graphene since the extremely small spin-orbit interaction (SOI) of carbon makes the 
spin dissipation that otherwise exists in most materials practically negligible [121 [13]. This 
suggests that whatever information is stored in the electronic spin of graphene can be re- 
tained for times considerably longer than those found in ordinary metals. Furthermore, this 
information can travel much longer distances with very little loss [TH4T6] . Not surprisingly, 
there is a growing interest in graphene-based spintronics as demonstrated by the volume of 
recent literature on the topic [17] . 

Driven by the necessity of a bandgap and by the growing interest in graphene-based 
spintronics, in this letter we propose a simple mechanism that not only produces a gapped 
electronic structure in graphene but that also spin-polarizes its current. We show that this 
effect arises quite simply by the combined presence of two key ingredients: the SOI and 
an externally applied magnetic field. While magnetic fields are controllable, the SOI of a 
material is normally constant and small in the case of carbon. Therefore, it might seem 
too ambitious to amplify both ingredients enough for the appearence of a possible gap. 
Nevertheless, recent discoveries have demonstrated that the SOI is enhanced when graphene 
is mechanically bent away from its planar geometry [T81 - I2T] suggesting that folding might 
function as a viable mechanism to induce a bandgap. In fact, here we show that folded 
graphene sheets in the presence of externally applied magnetic fields may display both 
a gapped electronic structure and spin-polarized currents. Not being possible to fold 3- 
dimensional structures, this possibility was never imagined with bulk materials but paves 



the way to a whole new approach of deahng with spin electronics in 2-diniensional systems, 
giving rise to the so-called origami spintronics. 

It is convenient to start by showing that the two aforementioned key ingredients indeed 
lead to a bandgap in graphene. Let us then consider a graphene sheet with SOI and in the 
presence of an externally applied magnetic field. For energies close to the Dirac point, the 
band Hamiltonian H is commonly described by a 2-dimensional Dirac equation. In this case, 
the Hamiltonian is given by [22]: H = vpl (g) a ■ p + 2\{ax ®Ty + ay® r^) + 7I (g) r^, where vp 
is the Fermi velocity, (J^^y^z are Pauli matrices acting on the sublattice space, t^. ^ ^ are Pauli 
matrices acting on the spin subspace, p = —iKV is the momentum operator, A represents 
the strength of the SOI and 7 is the Zeeman-field factor. Note that the Hamiltonian recovers 
the standard form for pristine graphene when both 7 and A vanish, so it is convenient to 
express the Hamiltonian H = Hq + V, where Hq is the pristine Hamiltonian and V accounts 
for the 7- and A-dependent perturbation. When expressed on the basis formed by the vector 
^ ~ ("^I' ^i' V'b ; V'ft )"^; where ip^ is the electron wavefunction on sublattice C with spin a, 
the Hamiltonian is written in matrix form as: 



li [Kx, Ky 



I —7 kx + iky —2iX \ 

7 kx + iky 

k,j. — iky —7 



(1) 



y 2i\ kx — iky 7 / 

The four corresponding eigenvalues are: 



E{k) = ±y P + ^2 + 2A2 ± 2^/\^ + k^{-f^ + A2) , (2) 

where k"^ = k^ + ky. Figure Hia) shows the dispersion E(k) for arbitrary values of A and 
7. Note the distinctive energy gap A = , ^'^ around E = 0. The characteristic linear 
dispersion relation of pristine graphene reappears from Eq. (|2| when 7 = A = 0. It is 
worth pointing out that either 7 or A alone is not capable to induce the bandgap opening, 
as shown in Figs. W[h) andfiyc). Only when both quantities are non-zero, will the bandgap 
appear. This is explained by the fact that 7 and A define two complementary energy scales 
and the magnitude of the bandgap A is determined ultimately by the smallest of these two 
quantities. 

Let us now consider folding the graphene sheet as a SOI-enhancing mechanism. In this 
case, however, the SOI is not homogeneously distributed but spatially hmited to the region 



of the graphene sheet surrounding the non-planar deformations. We consider a corrugated 
structure where curved deformations to the graphene sheet are ordered in a periodic struc- 
ture. This can be achieved by patterning a substrate in the form of periodic trenches [2S] 
on top of which a graphene sheet is deposited. In this way, the graphene sheet alternates 
between curved and flat regions, the dimensions of which will depend on the geometry of the 
patterned substrate. Figure [2]^a) shows an image of a typical corrugated sheet composed of 
a sequence of parallel ridges and troughs. While the diagram depicts a sheet of width L, it is 
supposed to represent a system where L ^ S, meaning that for practical purposes L — > oo. 
For the sake of simplicity, we consider the corrugated graphene sheet made of a series of 
curved regions, all composed of semi-circular cross sections with curvature radius R, equally 
spaced by a distance D, the cross section of which is schematically shown in Fig. |2)d. It is 
worth noting that this setup is equivalent to an array of half nanotubes joined seamlessly 
together by fiat nanoribbons. Since the material has to comply during folding, this kind of 
distortion of the graphene sheet does not introduce the disorder that is so deleterious to the 
transport properties in other gap-opening proposals. 

Differently from the previous case where we looked at how the electronic band structure 
changes with the inclusion of the SOI plus the Zeeman term across the entire system, in the 
corrugated case we investigate how the conductance of the material responds to the addition 
of exactly the same ingredients, now in a periodic arrangement as shown in Fig. |2]d. In this 
case both these ingredients are added to the pristine Hamiltonian only where the graphene 
sheet has a non-zero curvature. The advantage of studying the transport response of the 
system is that we are able not only to identify the opening of a bandgap, which appears 
as a region of zero conductance, but we can simultaneously investigate the spin dependence 
of the transport response when the conductance is finite. In this way, we can study the 
appearance of a bandgap and assess its potential impact on spintronics applications at the 
same time. 

Rather than using the Dirac equation to describe the electronic structure of the sys- 
tem, we consider the tight-binding (TB) Hamiltonian whose accuracy is not limited to 
a narrow energy range surrounding the Dirac point [26j. The Hamiltonian is written as 
H = X] ^i,o-^j',o- + V, where the operators Oj^^ io,j,a) creates (anihilates) an electron at site 
j and spin a on the sublattice A. The b, ^ and bj^^ are the corresponding operators for the 
B sublattice. The operator V = iX ^^^, X]a« ^ ' i^ap x diii)(yi^cii p — 7 Yliia ^^^a^ia represents 



the effects of a Rasliba spin-orbit coupling [22] plus a Zeeman field, dui is the position of 
site /' relative to site /. a is the vector formed by the Pauli matrices and the index a in the 
last sum takes the value +1 for majority spins and -1 for minority spins. The Rashba term 
arises whenever inversion symmetry in the direction perpendicular to the graphene sheet is 
broken, a condition realized for samples deposited on a substrate or in the presence of an 
electric field perpendicular to the sheet plane. 

The parameters A and 7 are dependent on the curvature and on the applied magnetic field, 
respectively. The former can be estimated from References fTIj [T5]. which studied the SOI 
in carbon nanotubes, and the latter can be extracted from characteristic values of magnetic 
fields. For the case shown in Fig. |2p the potential V acts only where the curved regions are 
located. In other words, we must consider an infinite long sheet of graphene with an array of 
M equally spaced regions of width ttR separated by a distance D. Following Refs. [121 US], 
A is inversely proportional to the radius R, which means that higher curvatures correspond 
to narrower curved sections with correspondingly stronger values for A. In the extreme limit 
of very high curvatures {R — )■ 0) the width of the regions with non-zero SOI becomes small, 
giving rise to a geometry depicted by Fig. ^. In all cases, we assume that the corrugated 
structure is connected by leads made of semi-infinite pristine graphene. 

The conductance is calculated by the Kubo formula expressed in terms of the single 
particle Green function G, which in turn, is defined as G{E) = {El — H)^^, where E is the 
energy and 1 is the identity operator. The zero bias DC conductance F is written as: 

F = ^ReTr ([G]ootoi[G]ntio - toi[G]iotoi[G]io) (3) 

Here, Gia^vp = [Gia^vp — (G'i//3,/Q,)*]/2z and [G]ab is the matrix formed by the Green function 
elements connecting unit cells A and B. a and /3 are spin indices. We define unit cells as 
lines along the "armchair" direction. The trace is taken over both site and spin indices. All 
Green functions in the above expression are evaluated at the Fermi level. 

In Fig. [3] we show the results for the spin-dependent conductance as a function of the 
Fermi level Ep for M = 10, 15, and 50. The first thing worth noticing is the opening of a 
gap as M increases, as shown in the inset of Fig. p[a). The magnitude of the gap increases 
with M and tends to a saturation value when M ^ 1. In addition, outside the gapped 
energy region, the conductance becomes spin dependent, i.e., the transport properties for 
the t and J, spin electrons become completely asymetric, which is a basic requirement for 



exploring the use of any material in spintronics. Such an asymmetry in the spin-dependent 
conductance means that these corrugated systems can function as spin polarizers for the 
electronic currents in graphene. The inset of Fig. [stb) shows the spin-polarization of the 
current, defined as Isp = (/-[- — J^) /(/-)- + /j,), as a function of Ep whereas the inset of Fig. |3](c) 
depicts Isp as a function of A. Note the sign change of Isp on either side of the band gap. 
This feature is seen in all our results. 

Having carried out the numerical evaluation of Eq. ([s]) within the TB model and ob- 
served the appearance of the bandgap as well as the spin asymmetry in the conductance, it 
is instructive to assess the robustness and generality of these findings by using alternative 
methods. With that in mind, we introduce a simplified approach that is far less computa- 
tionally involved and yet capable of capturing the contribution from the two key ingredients, 
namely the SOI and the Zeeman term. For the sake of simplicity, we investigate the limit 
i? — )■ 0, shown schematically in Fig. [2Fc). In this limit, the SOI-enhancing curvature becomes 
extremely large at the same time as the extension of the curved region becomes vanishingly 
small. The smooth ondulations shown in Figs. |2](a) and W[h) give way to equally spaced 
creases on the graphene sheet, ideally represented by Dirac delta functions that mark the 
locations where both the SOI and the Zeeman term are non-zero. Moreover, a further sim- 
plification is to assume electrons moving along the direction normal to the creases. The 
Hamiltonian in this case becomes H = IIq + X]7=i ^ ^{^ ~ JD)^ where the operators Hq and 
V have been previously defined, the only difference being that ky = and k = k^. In other 
words, the problem in question is now solved by considering Dirac equation for massless 
particles under the action of a perturbing potential made of a series of A- and 7-dependent 
delta functions. 

Care must be taken when solvig Dirac equation in the presence of delta function potentials 
because the wavefunctions become discontinous in this case [271 I2H] • The boundary condition 
needed to specify the wavefunctions on opposite sides of a given delta function at an arbitrary 
position, say a; = 0, is given by ^_|_ = e^^_, where A = —Vi^Ho)"^, \E'+ (^-) being the 
wavefunction on the right (left) hand side of the delta function position. So, for an array of 
M delta functions equally spaced by a distance D in between two leads, we write "^r = T'^l, 
where "^r (^l) is the wavefunction on the right (left) lead and T is the transfer matrix. 



The transfer matrix is written as T = [Z/(x/j)] ^ [e^ S{D)]^^ ^ L{xl), where: 



/^ikx g _^-ikx g \ 



Lix] 






^ikx 





^-ikx 
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^ikx 





^-ikx 



(4) 



and 



S{x) 



(5) 



/ cos kx i sin kx \ 
cos kx i sin /ex 
i sin A;x cos kx 
y i sin kx cos A;x y 
Xl and X/j are the positions marking the start of the left and right leads, respectively. 

It is now straightforward to extract the spin-resolved transmission coefficient from the 
transfer matrix T and plot it as a function of the energy, which is proportional to the 
wavevector k. Since the conductance is also proportional to the transmission coefficient, 
we can now compare the results of this simplified approach with the TB-based results. In 
Fig. |4] we show results for M = 2 delta functions a distance D apart calculated with the 
simplified approach as well as with the TB adapted for normal incidence. Quantitative 
agreement is not expected because the discrete delta function representation in the TB- 
model is fundamentally different from the continuous one. From a qualitative point of 
view, however, both methodologies lead to the appearance of a low-energy conductance gap 
followed by a spin-polarized current for energies outside the bandgap. This commonality 
found with two distinct methodologies is indicative of the robustness of our findings, which 
is further confirmed by recent calculations based on Density Functional Theory results also 
displaying bandgap features when the SOI is enhanced by curvature effects |29j . 

In addition to folding, adsorbants that deform the pristine planar structure of graphene 
sheets also function as SOI-enhancers. Therefore, it is possible to combine them both such 
that SOI-enhancing adsorbants might get attached to mechanically curved regions which 
themselves possess enhanced spin-orbit coupling. While we are not able to quantify the 
precise contributions coming from each one of the different mechanisms, curvature and 
dopants acting jointly together will certainly amplify the overall effect of spin-polarization 
of the charge currents as well as the bandgap in the electronic structure of the system. 



In summary, we have shown that by adding the Zeeman term and the SOI to the Hamil- 
tonian of pristine graphene, we force the opening of a bandgap in its electronic structure 
and induce its charge current to become spin polarized. Both additions to the Hamiltonian 
are experimentally feasible and are mimicked by an externally applied magnetic field and by 
deforming the graphene sheet out of its planar geometry through folding. With the increas- 
ing degree of experimental control on both the chemistry and the geometry of nanoscaled 
surfaces, it is possible that by engineering the shape and composition of graphene sheets 
we can create spin-polarized currents in an energy-gapped material. Such an origami-like 
technique is likely to bring about a whole new range of spintronic features in 2-dimensional 
systems not possible in 3-dimensional structures. 
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FIG. 1: Dispersion relation E{k) obtained from the four eingevalues of the Hamiltonian H. (a) In 
general the electronic structure displays a distinctive bandgap A, depicted by the vertical arrow. 
The energy gap disappears when (b) 7 / 0; A = or (c) 7 = 0; A / 
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FIG. 2: Schematic diagrams of the corrugated geometry adopted by our graphene sheets, (a) 
Characteristic shape of corrugated sheets of graphene of width L with a periodicity S. For practical 
purposes we assume L — t- oo. (b) Cross section of the corrugated graphene sheet used in our 
calculations, whereby sections of curved graphene with curvature radius R and width ttR are 
spaced by flat regions of width D. For distinction, these two different regions are separated by the 
horizontal dashed lines, (c) In the limit i? — )• 0, the corrugated sheet becomes a series of equally 
spaced sharp creases separated by flat sections of graphene. 
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FIG. 3: Spin-dependent conductances of a corrugated sheet made of M curved regions of radius 
R ~ 6A separated by flat sections of graphene of width D = 20. Except for the insets, black (red) 
lines represent t (i) spin-polarized conductances, (a) M = 10; 6 = 5 x 10^^; A = 5 x 10^^. The 
inset shows how the conductance gap scales with M. Note that it tends to a saturation value when 
M » 1. (b) M = 15; 5 = 5 X 10~^; A = 5 x 10^^. The inset depicts the spin polarization of 
the current Igp as a function of Ep. Only results outside the bandgap are shown, (c) M = 50; 
(5 = 5 X 10~^; A = 5 X 10^^. Shown in the inset is the dependence of Igp on the SOI parameter A. 
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FIG. 4: Electronic transport calculations on creased graphene sheets represented schematically by 
Fig. ^c). Spin-dependent conductance results for M = 2 delta functions for electrons constrained 
to move perpendicularly to the creases calculated based on the TB-model (main figure) and on the 
Continuous model (inset) for D = a, a being the lattice parameter of graphene. Black (red) lines 
represent f (i) spin-polarized conductances. 



13 



